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Introduction 


The  University  of  Rhode  Island  (URI)  together  with  the  Japan  Agency  for  Marine -Earth 
Science  and  Technology  (JAMSTEC),  Seoul  National  University  (SNU)  and  the  Korea 
Ocean  Research  &  Development  Institute  (KORDI)  undertook  a  joint  effort  to  study  the 
time  variability  of  the  Kuroshio  in  the  East  China  Sea  (ECS).  The  URI/JAMSTEC  group 
measured  the  offshelf  segment  of  the  flow  (between  the  550  m  isobath  of  the  continental 
shelf  and  the  Ryukyu  Islands),  principally  using  pressure-sensor  equipped  inverted  echo 
sounders  (PIES),  current  meter  equipped  PIES  (CPIES),  and  ship-mounted  acoustic 
Doppler  current  profilers  (ADCP).  The  SNU/KORDI  group  principally  measured  the 
onshelf  segment  of  the  flow  with  moored  ADCPs.  The  two  groups  made  hydrographic 
measurements  (with  XCTDs  and  CTDs,  respectively). 

The  following  Data  Report  is  presented  in  two  parts.  Part  I  focuses  on  the  data  collected 
by  the  URI/JAMSTEC  group  and  Part  II  focuses  on  the  data  collected  by  the 
SNU/KORDI  group. 
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Part  I.  The  URI/JAMSTEC  Measurements 


Introduction 

In  order  to  study  the  time  variability  of  Kuroshio  transport  and  structure  in  the  East  China 
Sea  (ECS),  two  lines  of  inverted  echo  sounders  (IES)  extending  across  the  axis  of  the 
Kuroshio  from  the  continental  shelf  to  the  Ryukyu  Island  chain  were  deployed  for  nearly 
two  years.  This  part  of  the  report  describes  the  data  recorded  by  these  instruments  and  its 
preliminary  processing.  During  the  deployment  cruise,  XCTD  casts  were  taken  at  each 
IES  site;  these  data  are  also  reported  along  with  ship  ADCP  data  taken  during  the  IES 
deployment  cruise.  In  addition,  the  technique  used  to  generate  GEM  fields  from  historical 
hydrographic  data  in  the  ECS  is  explained.  These  GEM  fields  will  be  used  to  infer  the 
dynamics  of  the  region  from  the  IES  acoustic-travel-time  (x)  records. 


The  ECS  study 

Instrumentation 

In  December  2002,  1 1  IESs  were  deployed  from  the  R/V  Yokosuka  of  the  Japan  Agency 
for  Marine-Earth  Science  and  Technology  (JAMSTEC).  In  November  2004,  all  were 
recovered  aboard  the  T/V  Kagoshima-maru  of  the  Faculty  of  Fisheries,  Kagoshima 
University.  Two  types  of  IES  instruments  were  deployed:  six  “CPIES”  instruments 
(manufactured  by  and  belonging  to  the  University  of  Rhode  Island)  and  five  “PIES” 
instruments  (manufactured  by  WHISL  and  belonging  to  the  Naval  Research  Laboratory). 
All  these  instruments  measured  round-trip  acoustic  travel  time  (t)  from  the  bottom  to  the 
surface,  together  with  near-bottom  pressure  and  temperature.  In  addition,  the  CPIESs 
were  equipped  with  Aanderaa  current  sensors  which  measured  current  velocity  and 
temperature  51m  above  the  sea  floor.  The  IESs  all  emitted  24  acoustic  pings  per  hour, 
from  which  a  single  hourly  acoustic  travel  time  value  was  calculated  (see  the  Data 
Processing  section  below),  but  the  timing  and  frequency  of  these  pings  were  different  in 
the  two  instrument  types  (see  Figure  1).  The  CPIESs  emitted  sets  of  4  12-kHz  pings 
every  10  minutes  over  an  entire  hour.  The  PIESs  sent  24  10-kHz  pings  at  10-second 
intervals  at  the  beginning  of  each  hour.  Both  the  CPIES  and  PIES  instruments  were 
equipped  with  Paroscientific  pressure  sensors  which  measured  hourly  near-bottom 
pressure  and  temperature  values.  In  the  CPIESs,  both  the  pressure  and  temperature  were 
measured  over  a  16  s  interval  following  the  first  burst  of  4  pings.  In  the  PIESs,  pressure 
was  averaged  over  the  entire  hour  and  temperature  was  averaged  over  the  last  minute  of 
the  hour. 

Each  CPIES  was  connected  through  a  50  m  cable  to  an  Aanderaa  Model  3820R  Doppler 
Current  Sensor,  which  recorded  current  speed  and  direction  as  well  as  temperature. 


Instrument  locations 

The  Kuroshio  flows  through  the  strait  east  of  Taiwan  into  the  ECS  and  along  the  edge  of 
the  continental  shelf  before  exiting  the  ECS  through  the  Tokara  Strait.  Our  IES  array, 
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north  of  Okinawa,  was  situated  in  the  Okinawa  Trough  between  the  shelf  and  the  Ryukyu 
Island  chain.  Instrument  locations  are  listed  in  Table  1  and  shown  in  Figure  2.  The  6 
CPIESs  were  deployed  across  the  axis  of  the  Kuroshio  from  the  ECS  continental  shelf 
slope  southeast  across  the  Okinawa  Trough  to  the  Ryukyu  Islands,  with  spacing  between 
instruments  ranging  from  ~20  to  ~40  km  (instruments  were  spaced  closest  beneath  the 
Kuroshio).  The  CPIESs  were  placed  near  the  PN  line,  along  which  hydrographic  data  are 
regularly  collected  four  times  per  year  by  the  Nagasaki  Marine  Observatory,  Japan 
Meteorological  Agency.  The  5  PIESs  were  deployed  across  the  Kuroshio  axis  ~40  km 
downstream  of  the  PN  line  (towards  the  northeast).  The  PIES  spacing  was  also  between 
~20  and  ~40  km. 


CPIES 


tau 

press. 

temp. 

Aanderaa  current  &  temp. 


PIES 


tau 

press. 

temp. 
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20  30  40 

time  (minutes) 


50 
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Figure  1.  Timing  of  measurements  in  CPIES  and  PIES  instruments. 


Table  1.  IES  instrument  depths  and  locations.  Spacing  is  separation  in  x-direction  along  128  °T. 


Instrument 

C  =  CPIES 

P  =  PIES 

IES 

Serial 

No. 

Current 
Sensor  Serial 
No. 

Latitude 

(°N) 

Longitude 

(°E) 

Spacing 

(km) 

Pressure 

level 

(dbar) 

Cl 

36 

164 

28.3425 

127.0058 

19.43 

19.07 

31.81 

36.35 

41.85 

553.5 

C2 

43 

165 

28.2622 

127.1822 

1005.7 

C3 

45 

166 

28.1234 

127.2968 

1131.4 

C4 

53 

168 

27.9706 

127.5711 

1338.9 

C5 

63 

169 

27.7712 

127.8644 

1115.5 

C6 

68 

170 

27.5369 

128.1973 

891.2 

PI 

117 

n/a 

28.6370 

127.2637 

18.95 

20.25 

32.71 

37.90 

1033.1 

P2 

118 

n/a 

28.5273 

127.4124 

1121.3 

P3 

106 

n/a 

28.3533 

127.6808 

856.3 

P4 

111 

n/a 

28.1743 

127.9462 

1079.9 

P5 

116 

n/a 

27.9614 

128.2481 

749.9 

Data 

With  the  following  exceptions,  acoustic-travel-time,  pressure,  temperature  and  current 
data  records  are  all  of  excellent  quality  and  complete.  (1)  Current  data  from  the 
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Aanderaa  sensor  on  instrument  C6  are  intermittent  from  April  2004  until  the  end  of  the 
record.  It  is  not  clear  what  caused  this  failure.  (2)  After  about  one  week,  the  current 
sensor  on  C5  failed.  Also  about  1/3  is  missing  of  the  Paroscientific  pressure  and 
temperature  records  from  this  instrument  (though  the  acoustic-travel-time  record  is 
complete).  These  problems  on  C5  appear  to  have  been  caused  by  a  faulty  o-ring  seal  on 
the  current-sensor  connector. 


Figure  2.  Map  of  instrument  locations.  CPIES  sites  are  designated  C;  PIES  sites  are  designated  P;  bottom- 
mounted  ADCP  sites  (see  Part  II  of  this  report)  are  designated  A. 


Data  Processing 

Removing  the  “jumps” 

The  pressure  records  of  the  two  shallowest  instruments  show  several  pressure  “jumps,” 
probably  because  the  instruments  were  dragged  by  bottom-fishing  gear.  The 
westernmost  CPIES  (Cl  at  about  550  m  depth)  was  moved  twice,  leaving  the  instrument 
about  1  m  deeper  after  the  first  incident  and  another  5  m  deeper  after  the  second.  The 
easternmost  PIES  (P5  at  760  m  depth)  was  moved  once  about  1  m  deeper. 

Before  processing  the  data  from  these  two  instruments,  their  pressure  and  acoustic-travel- 
time  records  were  modified  to  remove  the  effects  of  these  changes  in  instrument  depth. 
Each  jump  was  corrected  by  adding  to  the  entire  record  prior  to  the  jump  an  amount  equal 
to  the  size  of  the  jump. 

A  preliminary  detiding  of  the  data  records  containing  jumps  was  first  performed  using  the 
response  method  (Munk  and  Cartwright,  1966).  The  time  of  the  pressure  jumps  was 
identified  from  this  preliminary  detided  record.  The  magnitude  of  each  jump  was 
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determined  by  differencing  the  mean  detided  pressures  24  hours  before  and  after  the 
jump  (Table  2).  Once  the  magnitudes  of  the  pressure  jumps  were  detennined  and 
converted  to  depth  changes,  the  magnitudes  of  the  jumps  in  the  acoustic-travel-time 
records  were  calculated  using  the  average  sound  speeds  for  the  depth,  temperature  and 
salinity  of  each  instrument  (the  temperature  and  salinity  were  determined  from  the  mean 
of  ECS  hydrocast  data  at  the  nominal  IES  depths).  The  original  (i.e.,  non-detided) 
pressure  and  acoustic-travel-time  records  were  then  corrected  for  these  jumps  (see  Table 
2).  Typical  vertical  temperature  gradients  at  Cl  and  P5  are  0.01°C/m,  so  the  sudden 
increases  in  pressure  would  have  resulted  in  temperature  jumps  of  only  a  few  hundredths 
of  a  degree.  We  have  not  made  any  temperature-jump  corrections. 


Table  2.  Jump  corrections.  Note:  yearday  is  in  days  since  2002  January’  1,  0  GMT. 


Instrument 

Beginning  of  jump  (yearday) 

AP  (dbar) 

At  (s) 

Cl 

505.1671 

0.7 

0.0009 

Cl 

860.8338 

4.7 

0.0063 

P5 

621.1489 

0.6 

0.0008 

Pressure  and  temperature 

Following  jump-removal,  each  pressure  record  was  again  detided  using  the  response 
method.  While  the  precision  of  the  pressure  transducers  is  good,  they  tend  to  drift  over 
time,  particularly  when  stressed  at  high  pressures.  The  instruments  in  the  ECS  study 
were  fairly  shallow  (<1500  m)  and  thus  have  relatively  small  drifts.  Drifts  were 
determined  by  fitting  exponential-plus-linear  curves  (see  Table  3)  to  the  dejumped  and 
detided  pressure  records.  These  drifts,  shown  in  Figure  3,  were  then  subtracted  from  the 
pressure  records.  The  resulting  detided,  dedrifted  records  were  lowpass  filtered  using  a 
2nd  order  Butterworth  filter,  run  forward  and  backward,  with  cutoff  period  of  48  hours, 
and  were  subsampled  at  12-hour  intervals  (Figure  4).  The  Paroscientific  temperature 
records  were  despiked,  and  then  filtered  and  subsampled  in  the  same  manner  (Figure  5). 
Comparison  of  temperatures  measured  by  the  Paroscientific  sensors  to  hydrographic  data 
and  to  temperatures  measured  by  the  Aanderaa  current  sensors  indicate  that  the 
Paroscientific  temperature  measurements  have  temperature  offsets.  For  PIES 
instruments,  temperature  measured  by  the  Paroscientific  gages  is  approximately  0.5°C 
higher  than  the  temperature  one  would  predict  from  historical  hydrocasts.  For  CPIES  this 
offset  varies  greatly  between  instruments  (from  0.1 5°C  colder  than  hydrocast  data  to 
1.6°C  warmer).  Thus  only  temperature  variation,  rather  than  absolute  temperature,  is 
reliably  measured  by  the  Paroscientific  gages  on  the  CPIESs  and  PIESs. 

Acoustic  travel  time 

As  described  in  the  Instrumentation  subsection  above,  the  IES  recorded  24  acoustic  travel 
times  each  hour.  These  were  converted  to  a  single  hourly  acoustic-travel-time 
measurement  using  the  “quartile  method”  described  as  follows.  For  a  given  hour  h,  the  t 
samples  are  first  passed  through  two  stages  of  windowing  in  which  samples  are 
eliminated  if  far  removed  from  either  the  median  of  the  medians  of  the  neighboring  hours 
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or  from  the  first-quartile  point  of  the  hour-/?  samples.  From  the  remaining  acceptable 
n< 24  samples,  the  x  value  assigned  to  hour-/?  is  taken  to  be  the  average  of  the  ??/6  samples 
beginning  with  the  last  sample  of  the  first  quartile.  For  instance,  for  hour  /?,  if  ??=24  (i.e., 
all  samples  in  the  hour  are  acceptable),  the  average  of  the  6th,  7th,  8th,  and  9th  smallest  x 
samples  is  assigned  as  the  hour-/?  x  value. 

The  resulting  records  were  despiked,  then  lowpass  filtered  using  a  2nd  order  Butterworth 
filter,  run  forward  and  backward,  with  cutoff  period  of  48  hours,  and  finally  subsampled 
at  12  hour  samples.  These  processed  acoustic-travel-time  series  are  plotted  in  Figure  6. 


Table  3.  Coefficients  for  the  pressure  sensor  drift  equation:  drift  (dbar)  =  AeBt  +  Ct  +  D,  t  is  time  in  days 
since  the  beginning  of  the  record.  Note  yearday  is  in  days  since  2002  January >  1,  0  GMT. 


IES 

t  =  0  time 
(yearday) 

A 

B 

C 

D 

Cl 

339.3754 

-0.13132 

-0.04085 

0.00030 

553.51 

C2 

339.3341 

0 

0 

0.00002 

1005.75 

C3 

339.5430 

-0.02373 

-0.04945 

0.00002 

1131.43 

C4 

339.6258 

-0.04468 

-0.07860 

0.00000 

1338.86 

C5 

339.7090 

-0.05135 

-0.06790 

0.00000 

1115.48 

C6 

339.9173 

-0.17442 
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0.00004 
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PI 
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0.00007 
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-0.02220 
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P5 
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0 

0 
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Figure  3.  Pressure  Sensor  Drifts.  Detided  pressure  data  (jumps  removed)  for  CPIESs  and  PIESs  shown  in 
blue.  The  computed  drifts  for  each  instrument  are  shown  in  red. 
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CPIES 


Figure  4.  CPIES  and  PIES  pressure  time  series  (jumps  removed).  Records  have  been  detided,  dedrifted 
and  lowpass  filtered.  Means  have  been  removed  and  series  have  been  displaced  for  clarity  by  0.2  dbar 
(except  P5  which  has  been  displaced  by  0.4  dbar). 
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CPIES 


1/1/03  4/1/03  7/1/03  10/1/03  1/1/04  4/1/04  7/1/04  10/1/04 


Figure  5.  CPIES  and  PIES  temperature  time  series  (blue  lines).  Aanderaa  Current  Sensor  temperature 
time  series  (red  lines).  Records  have  been  despiked  and  lowpass  filtered.  Means  have  been  removed  and 

successive  series  displaced  for  clarity  by  1  °C. 
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-0.03 


1/1/03  4/1/03  7/1/03  10/1/03  1/1/04  4/1/04  7/1/04  10/1/04 

Figure  6.  CPIES  and  PIES  acoustic-travel-time  time  series  (jumps  removed).  Records  have  been  despiked 
and  lowpass  filtered.  Means  have  been  removed  and  successive  series  displaced  for  clarity  by  5  ms. 
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Current  data 

Current  sensor  records  were  corrected  for  the  local  magnetic  declination  at  each  site  (data 
from  http://www.ngdc.noaa.gov/seg/geomag/declination.shtml;  roughly  5°W),  so  U  is 
eastward  and  V  is  northward.  The  velocities  were  then  lowpass  filtered  using  a  2nd  order 
Butterworth  filter,  run  forward  and  backward,  with  cutoff  period  of  48  hours  The 
resulting  records  were  subsampled  at  12-hour  intervals  and  are  shown  as  stick-plots  in 
Figure  7.  The  partial  record  from  the  current  sensor  on  C5,  which  failed  after  6  days,  is 
shown  in  Figure  8.  Table  4  gives  the  basic  statistics  for  the  components  of  the  current 
vectors  ahown  in  Figures  8  and  9  in  the  directions  along  the  array  (128°T)  and  across  the 
array  (38°T).  The  temperatures  recorded  by  the  Aanderaa  Current  Sensors  are  shown  in 
Figure  5  together  with  the  respective  CPIES  temperature  records,  except  for  site  C5, 
where  the  Aanderaa  Current  Sensor  failed. 


20_i_ i_ i_ i_ i_ i_ i_ i 
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Figure  7.  Current  time  series.  Vertical  axis  is  toward  True  north. 


Figure  8.  Partial  current  record  for  C5.  Vertical  axis  is  towards  True  North. 
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Table  4.  Basic  Statistics  for  CPIES  measurement  of  currents  51  m  above  bottom  for  Ur  (towards  128  °)  and 
Vr  (towards  38  °)  components  (48-hour  lowpass  fdtered).  Max  and  Min  denote  maximum  and  minimum 
values  of  current  components.  The  direction  of  vector  current  is  measured  clockwise  from  True  North. 


CPIES 

Depth 

(m) 

Vel. 

Mean 

(cm/s) 

STD 

(cm/s) 

Max. 

(cm/s) 

Min. 

(cm/s) 

Vector  Mean 

Data 

Returned 

(%) 

Dir. 

(°) 

Speed 

(cm/s) 

Cl 

498 

Ur 

1.42 

3.51 

18.09 

-10.49 

172.4 

1.99 

100 

Vr 

-1.39 

7.83 

26.07 

-34.00 

C2 

946 

Ur 

0.46 

3.14 

12.82 

-9.57 

208.9 

2.92 

100 

Vr 

-2.88 

5.56 

16.13 

-21.18 

C3 

1070 

Ur 

-4.33 

3.44 

7.50 

-20.55 

16.6 

11.88 

100 

Vr 

11.06 

6.80 

29.06 

-8.15 

C4 

1275 

Ur 

1.29 

2.71 

21.12 

-7.06 

146.3 

1.36 

98.6 

Vr 

-0.43 

2.84 

15.28 

-9.01 

C5 

1054 

Ur 

-1.26 

0.47 

-0.81 

-2.25 

241.1 

3.21 

1.0 

Vr 

-2.95 

0.45 

-2.21 

-3.74 

C6 

832 

Ur 

-1.41 

3.16 

10.67 

-20.65 

236.2 

4.51 

85.6 

Vr 

-4.28 

2.92 

5.67 

-19.37 

Deployment  cruise  data 

XCTD  data 

XCTD  casts  were  taken  after  each  instrument  deployment  during  the  December  2002 
cruise  aboard  the  R/V  Yokosuka.  Figure  9  shows  the  temperature  and  salinity  sections  for 
each  of  the  two  lines  of  instruments.  Table  5  lists  the  locations  and  depths  of  these  casts. 

Near-surface  current  data 

Ship  speed  relative  to  sea  water  measured  by  electromagnetic  log  and  ship  heading 
measured  by  gyro-compass  were  collected  every  5  seconds  during  the  deployment  cruise 
aboard  the  R/V  Yokosuka  together  with  water  temperature  at  the  ship’s  hull  bottom  (6  m 
depth)  and  ship’s  DGPS  position.  The  30-minute-mean  near-surface  current  vectors 
calculated  from  these  data  are  plotted  in  Figure  10  together  with  the  near-surface 
temperature.  Due  to  characteristics  of  the  above-mentioned  conventional  measurement 
system,  the  calculated  current  vectors  may  have  large  errors  during  times  when  the  ship 
was  not  running  at  high  speed. 
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Figure  9.  Hydrographic  sections  across  the  CPIES  and  PIES  lines.  Open  circles  show  the  longitude  and 

depth  of  nearby  IES  instruments 


.  Table  5.  Locations,  depths  and  times  of  XCTD  casts  during  2002  deployment  cruise. 


Nearest 

IES 

Longitude 

(°E) 

Latitude 

(°N) 

Cast 

depth  (m) 

Cast  time 
mo/day  hr  GMT 

Cl 

127.0065 

28.3473 

518 

12/5 

21 

C2 

127.1794 

28.2599 

990 

12/5 

23 

C3 

127.2989 

28.1217 

741 

12/6 

0 

C4 

127.5752 

27.9699 

1001 

12/6 

2 

C5 

127.8661 

27.7712 

1001 

12/6 

4 

C6 

128.1997 

27.5361 

886 

12/6 

6 

PI 

127.2619 

28.6367 

1001 

12/6 

21 

P2 

127.4120 

28.5255 

1001 

12/6 

23 

P3 

127.6785 

28.3526 

858 

12/7 

1 

P4 

127.9439 

28.1740 

1001 

12/7 

3 

P5 

128.2464 

27.9610 

719 

12/7 

5 
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SST  and  SFC  on  Dec. 5-7,  2002 
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Figure  10.  Near-surface  currents  and  sea  surface  temperatures  (SST)  during  the  deployment  cruise. 
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Recovery  cruise  data 

Ship  ADCP  data 

The  current  velocities  at  50  depths  beginning  at  26  m  with  16  m  intervals  were  collected 
every  second  by  the  75  kHz  Ocean  Surveyor  Vessel-Mounted  ADCP  (RDI)  during  the 
recovery  cruise  aboard  the  T/V  Kagoshima-maru,  together  with  water  temperature  at  the 
ship’s  hull  bottom  (5  m  depth)  and  ship’s  DGPS  position.  The  misalignment  angle  of  the 
ADCP  system  was  detennined  to  be  0.8996  degree  anti-clockwise  from  the  available 
bottom  tracked  data  during  the  cruise.  From  the  short  term  (1  minute)  average  (STA)  of 
ADCP  data  and  navigation  data,  the  10-minute  mean  current  vectors  at  26  m  depth  were 
calculated  and  are  plotted  in  Figure  1 1  together  with  water  temperature  at  5  m  depth. 


126.5  127  127.5  128  128.5  129  129.5  130 

Longitude  (E) 


Figure  11.  Near-surface  currents  and  sea  surface  temperatures  (SST)  during  the  recovery’  cruise. 


GEM 

The  vertical  structures  of  the  current  and  temperature  fields  are  related  to  vertically 
integrated  quantities  such  as  acoustic  travel  time  (e.g.,  a  deeper  thennocline  is  associated 
with  a  shorter  acoustic  travel  time).  For  a  given  region,  this  relationship  (if  it  exists) 
between  vertical  temperature  profile  (or  specific-volume-anomaly  profile)  and  acoustic 
travel  time  can  be  determined  from  historical  hydrographic  data.  That  part  of  the  vertical 
structure  captured  by  this  relationship  is  the  Gravest  Empirical  Mode  (GEM). 
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To  compute  GEMs  for  the  region  we  used  1833  historic  ECS  hydrographic  profdes 
obtained  from  North  Pacific  Hydrobase  (Macdonald  et  ah,  2001)  and  the  Nagasaki 
Marine  Observatory  (NMO),  Japan  Meteorological  Agency.  All  of  these  profiles 
extended  to  at  least  700  dbar.  Data  were  extracted  at  10  dbar  intervals.  The  data  were 
quality  controlled  by  removing  those  hydrocasts  having  density  inversions.  For  this 
purpose,  if  the  density  at  any  pressure  level  was  0.04  kg/m  '  lower  than  the  density  in  the 
overlying  layer  (i.e.,  10  dbar  shallower),  it  was  identified  as  a  density  inversion. 

In  many  applications,  a  single  GEM  field  for  a  particular  parameter,  such  as  temperature, 
is  used  as  a  lookup  table  for  a  given  region.  In  the  ECS,  however,  it  was  necessary  to 
construct  two  separate  fields:  a  localized  GEM  for  a  small  region  around  C6  and  a  main 
GEM  for  the  rest  of  the  IES  sites.  Detailed  bathymetry  data  (Choi  et  ah,  2002)  show  that 
C6  (885  m  depth)  was  isolated  from  the  rest  of  the  ECS  because  it  was  situated  within  a 
small  Y-shaped  basin  with  a  sill  depth  of  -800  m  and  a  maximum  bottom  depth  of 
-1200  m  (Figure  12).  Being  just  off  the  west  coast  of  the  island  Okinoerabushima  in  the 
Ryukyu  Island  chain,  we  call  this  the  Okinoerabushima  Basin.  A  comparison  of 
temperature  profiles  from  the  1833  hydrographic  stations  shows  that  water  within  this 
basin  is  warmer  than  water  at  the  same  depths  elsewhere  in  the  ECS,  typically  about 
1.5°C  warmer  at  850  dbar.  The  basin  appears  to  be  completely  filled  with  Pacific 
Intermediate  Water.  Only  the  1720  ECS  profiles  from  outside  the  basin  are  used  to 
compute  the  ECS  main  GEM,  while  the  localized  GEM  is  calculated  exclusively  from  the 
113  “warm”  hydrocasts  taken  in  the  basin. 


1 27°E  128°E 


|  I  I  I  I  |  I  I  I  I  |  I  I  I  I  |  I  I  I  I  |  I  I 

0  0.5  1  1.5  2 

Depth  (km) 

Figure  12.  Detailed  map  of  the  region  around  C6.  Blue  contour  indicates  the  800  m  isobath.  Dots  show 
IES  positions.  Bathymetry’  data  from  Choi  et  al.,  2002. 
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To  compute  the  main  and  localized  GEMs  the  hydrographic  data  for  each  of  the  two 
regions  were  first  de-seasoned  using  the  SLACTS*  method  of  Book  (1998),  with  the  year 
split  into  5  bins  and  the  center  of  the  bins  shifted  8  times.  Figure  13  shows  the  SLACTS 
curves  for  temperature  and  specific-volume-anomaly  signals  at  different  depths.  Vertical 
round-trip  acoustic  travel  time  from  700  dbar  to  the  surface,  x7oo,  was  calculated  from  the 
historical  hydrographic  data  using  a  constant  value  for  gravitational  acceleration,  g  = 
9.8  m/s2  (rather  than  the  local  g).  This  was  then  deseasoned  using  the  SLACTS  method; 
in  this  case  (since  only  casts  beyond  700  dbar  could  be  used)  the  year  was  split  into  4 
bins  with  the  center  of  the  bins  shifted  10  times.  The  x70o  SLACTS  curve  is  shown  in 
Figure  14.  A  spline  fit  was  used  to  quantify  the  relationship  between  the  deseasoned 
hydrographic  temperature  data  and  deseasoned  hydrographic  x7oo  at  each  pressure  level 
from  the  surface  to  2000  dbar  in  10  dbar  increments.  From  this,  we  generated  a  lookup 
table  of  temperature  as  a  function  of  pressure  (in  10  dbar  intervals)  and  x7oo  (in  0.1  ms 
intervals).  The  temperature  and  specific-volume-anomaly  GEMs  computed  from  these 
deseasoned  data  are  shown  for  the  main  and  localized  region  in  Figure  15  and  Figure  16, 
respectively. 

Seasonal  Signals 

’cold’  casts  ’warm’  casts 


Figure  13.  SLACTS  Temperature  (top)  and  Specific-Volume- Anomaly  (bottom)  seasonal  signals,  from  the 
surface  (largest  amplitude)  to  130  dbar  in  10  dbar  increments  for  the  ‘cold’  casts  (left)  and  ‘warm  ’  casts 

(right). 


“SLACTS”  stands  for  “Surface  Layer  Annual  Correction  for  Temperature  and  Salinity”. 
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'cold'  casts  ’warm’  casts 


Figure  14.  SLACTS  acoustic-travel-time  seasonal  signal  for  ‘cold’  casts  (left)  and  ’warm  ’  casts  (right). 


In  order  to  evaluate  how  well  the  GEM  represents  the  vertical  structure  of  the  water 
column,  the  GEM-predicted  temperature  and  specific-volume-anomaly  profiles  are 
compared  to  the  actual  temperature  and  specific-volume-anomaly  profiles  of  the 
historical  hydrographic  data.  For  each  hydrocast,  X700  is  calculated  from  the  data.  Then 
this  is  used  to  look  up  the  GEM-predicted  temperature  and  specific -volume-anomaly 
profiles.  The  rms  error  between  the  GEM-predicted  and  hydrocast  values  is  calculated. 
Then  the  average  rms  error  in  0.5  ms  by  10  dbar  bins  is  calculated  and  contoured  as  a 
function  of  pressure  and  X700. 
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Figure  15.  Main  ECS  OEMs  and  error  fields.  Upper  panels  are  for  temperature,  lower  panels  are  for 

specific  volume  anomaly. 
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Figure  16.  Localized  ECS  OEMs  and  error  fields.  Upper  panels  are  for  temperature,  lower  panels  are  for 

specific  volume  anomaly. 
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In  order  to  use  the  GEM  lookup  table  which  is  based  on  acoustic  travel  times  referenced 
to  700  dbar,  the  travel  time  data  from  each  IES  instrument  must  be  converted  from  travel 
time  referenced  to  the  instrument’s  pressure  level,  tp,  to  travel  time  referenced  to 
700  dbar,  T700.  First  the  pressure  level  of  each  instrument  is  determined  from  the  mean  of 
its  dedrifted  pressure  record  minus  atmospheric  pressure  (10.1325  dbar  =  1  atm).  Then 
hydrographic  data  are  used  to  calculate  t7oo  and  xp  for  each  historic  ECS  hydrographic 
profile  extending  to  this  pressure  level.  These  are  plotted  one  against  the  other  in  Figure 
17.  Also  shown  in  this  figure  are  least-squares  fitted  straight  lines,  except  for  Cl  where  a 
second  degree  polynomial  is  fitted  since  it  provides  a  significantly  better  fit.  Table  6 
gives  the  coefficients  for  all  these  fitted  functions. 

Only  the  113  hydrocasts  taken  in  the  Okinoerabushima  Basin  were  used  to  determine  the 
conversion  factor  for  C6.  The  1720  other  ECS  hydrocasts  were  used  to  determine  the 
conversion  relationships  for  all  the  other  instrument  sites. 


xp  (ms) 


Figure  1 7.  x700  plotted  with  Tp  for  each  of  the  IES  sites.  Each  of  the  best-fit  lines  (red)  is  a  linear  fit, 
except  for  Cl  which  is  fitted  with  a  second-degree  polynomial  (Table  5). 
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Table  6.  Coefficients  for  the  conversion:  Tyoo=  cia  +  a  if  +  a2tp2  (r  in  ms),  representing  the  red  lines  in 

Figure  1 7. 


IES 

3o 

ai 

a2 

rms  error  (ms) 

Cl 

-12648.00 

36.78 

-0.024897 

0.30 

C2 

-274.40 

0.91 

0 

0.19 

C3 

-415.67 

0.90 

0 

0.23 

C4 

-638.96 

0.89 

0 

0.24 

C5 

-396.85 

0.90 

0 

0.22 

C6 

-163.97 

0.93 

0 

0.14 

PI 

-300.81 

0.91 

0 

0.20 

P2 

-405.32 

0.90 

0 

0.22 

P3 

-125.92 

0.94 

0 

0.12 

P4 

-353.30 

0.90 

0 

0.21 

P5 

-27.82 

0.97 

0 

0.05 

25 


Part  II.  The  SNU/KORDI  Measurements 


Introduction 

Two  bottom-mounted  Acoustic  Doppler  Current  Profiler  (ADCP)  moorings  were 
deployed  in  the  East  China  Sea  (ECS)  for  13  months  to  observe  the  part  of  the  Kuroshio 
flowing  over  the  outer  continental  shelf  of  the  ECS.  The  main  objectives  of  this  study  are 
to  quantify  the  entire  Kuroshio  transport  and  to  investigate  the  dynamics  of  its  temporal 
variability  using  the  ADCP  data  together  with  CPIES  (current-and-pressure-sensor- 
equipped  inverted  echo  sounder)  data  simultaneously  obtained  in  the  Okinawa  Trough. 
This  report  describes  the  data  recorded  by  the  ADCPs  and  the  preliminary  processing. 
During  the  4  ADCP  mooring  cruises,  CTD  casts  were  conducted  along  the  CPIES-ADCP 
mooring  line  and  these  data  are  also  reported. 


ADCP  moorings 


Instrumentation 

The  ADCPs  were  bottom  moored  at  two  sites  located  on  the  outer  continental  shelf  of  the 
ECS  (Figure  18,  Table  7).  The  ADCP  moored  at  the  shallower  site  (station  Al)  was 
housed  in  a  bamacle-shaped  trawl-resistant  bottom  mount  shown  in  Figure  19a  (Barny 
TRBM;  Perkins  et  al.,  2000).  The  TRBM  has  proven  to  be  effective  for  measuring 
current  profiles  in  active  fishing  areas  for  relatively  long  times,  e.g.,  up  to  6  months  in  the 
Korea/Tsushima  Strait  (Teague  et  al.,  2002).  Currently,  the  TRBM-ADCP  package  can 
only  be  used  in  areas  shallower  than  150  m  due  to  the  limited  length  of  the  recovery  line. 
Thus,  the  ADCP  mooring  at  the  deeper  site  (station  A2)  was  made  using  a  subsurface 
foam  buoy  as  shown  in  Figure  19b.  The  field  work  consisted  of  four  cruises  to  improve 
the  likelihood  of  sufficient  data  return,  despite  the  risk  of  instrument  loss  in  this 
intensively  fished  region:  initial  deployment,  two  turnarounds,  and  final  recovery  (Table 
8).  In  October  2003,  two  ADCP  moorings  were  initially  deployed  (Leg  I)  from  the  R/V 
Onnuri  of  the  Korea  Ocean  Research  and  Development  Institute.  Turnarounds  of  the 
moorings  (beginnings  of  Leg  II  and  III,  respectively)  were  made  in  November  2003  from 
T/V  Kagoshima-maru  of  the  Faculty  of  Fisheries,  Kagoshima  University,  and  in  May 
2004  from  R/V  Tamyang  of  Pukyung  National  University.  The  ADCP  moorings  were 
finally  recovered  aboard  R/V  Tamyang  in  November  2004. 

The  ADCPs  used  are  RDI  Model  WHS300,  utilizing  “Janus”  beam  geometry  (four 
beams)  with  300  kHz  frequency.  The  location,  mooring  depth,  and  setting  of  the  ADCPs 
for  each  leg  are  given  in  Table  7.  The  ADCPs  recorded  current  speed  and  direction  at 
elevations  from  5-10  m  to  about  120  m  above  the  sea  floor  at  4  or  8  m  intervals.  The 
ADCPs  also  recorded  pressure  and  temperature  at  their  mooring  depths.  The  ADCP 
mooring  depths  in  Table  7,  152m  at  Al  and  285  m  at  Al,  were  estimated  using  the 
measured  pressure  data  from  Legs  III  and  II,  respectively.  The  ADCPs  at  stations  Al  and 
A2  were  respectively  placed  at  about  0.5m  and  5  m  above  the  sea  floor.  Hence  the  total 
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water  depths  at  stations  A1  and  A2  are  about  152.5  m  and  290  m,  respectively.  As  the 
measurement  range  of  the  WH300  ADCP  is  only  about  120  m,  vertical  profiles  of 
currents  are  measured  in  only  the  deeper  parts  of  the  water  column,  about  80%  and  40% 
of  total  water  depths  at  stations  A1  and  A2,  respectively.  An  attempt  was  made  to 
measure  the  upper  currents  near  station  A2  by  deploying  an  additional  mooring  equipped 
with  an  ADCP  at  140  m  and  an  RCM-1 1  current  meter  at  165  m  during  Leg  III,  but  this 
mooring  was  lost. 

Correlation  and  error  velocity  were  set  to  their  recommended  default  values,  64  and 
2  m/s,  respectively,  for  all  ADCPs.  When  an  ADCP  is  collecting  data,  the  correlation  of 
each  beam  is  handled  independently  for  the  purpose  of  data  quality  analysis.  Correlation 
is  essentially  a  measure  of  how  much  the  particle  distribution  has  changed  between  phase 
measurements.  The  less  the  distribution  has  changed,  the  higher  the  correlation,  and  the 
more  precise  the  velocity  measurement.  The  ADCPs  used  here  record  the  velocity  as  bad 
if  the  individual  correlation  from  at  least  2  beams  is  below  a  pre-determined  threshold 
(64).  If  the  error  velocity  reported  is  greater  than  the  error  velocity  threshold  (2  m/s),  then 
the  data  is  also  marked  as  bad.  The  low  correlation  and  large  error-velocity  setting  are  to 
avoid  data  loss. 


Table  7.  Mooring  locations  and  depths,  and  settings  of  ADCPs  for  each  leg. 


ADCP 

A1  Leg  III 

A2  Leg  I 

A2  Leg  H 

A2  Leg  m 

Period  (GMT) 

07:00  14  May.  2004  - 

17:00  1  Oct.  2003  - 

04:30  8  Nov..  2003  - 

10:30  14  May.  2004  - 

23:00  3  Nov..  2004 

03:00  7  Nov..  2003 

07:00  14  May,  2004 

23:00  3  Nov..  2004 

Latitude  f"W) 

28.5577 

28.4072 

28.4057 

28.4063 

Longitude  (°E) 

126.6699 

126.9113 

126.9098 

126.9108 

Bin  length  (m) 

4 

8 

8 

4 

Water  depth  (m) 

152.5 

290 

290 

290 

ADCP  Mooring  depth  (m) 

152 

285 

285 

285 

Fir  st  bin  depth  (m) 

146 

275 

275 

279 

Pings  per  ensemble 

60 

180 

50 

60 

Time  per  ping  (sec) 

30 

10 

36 

30 

Ensemble  interval  (min) 

30 

30 

30 

30 

ADCP  Serial  No. 

2279 

2937 

3569 

2937 

Table  S.  ADCP  mooring  time  coverage. 


ADCP  Moorings 

Period 

Oct.,  03 

Nov. 

Dec. 

Jan.,  04 

Feb. 

Mar. 

Apr. 

May 

Jtm. 

Jul. 

Allg. 

Sep. 

Oct. 

Nov. 

A1  Leg  in 

A2  Leg  I 

A2  Leg  n 

A2  Le  g  El 

Instrument  locations 

To  observe  the  part  of  the  Kuroshio  over  the  outer  continental  shelf  of  the  ECS,  the 
ADCPs  were  bottom  moored  at  depths  shallower  than  300m.  The  ADCP  moorings  were 
placed  northwest  of  a  CPIES  array  near  the  PN  line,  along  which  hydrographic  data  are 
regularly  collected  by  the  Nagasaki  Marine  Observatory,  Japan  Meteorological  Agency. 
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0  (m) 


Figure  18.  Positions  ofADCP  moorings  (Al,  A2)  and  CTD  stations  in  the  East  China  Sea.  CTD  stations 
C1~C6  nearly  coincide  with  locations  of  CPIES  moorings.  Bathymetry’  is  based  on  ET0P02  in  meter  units 
(http://www.ngdc.noaa.gov/mgg/fliers/01mgg04.html). 


The  three  deployment  locations  at  A2  were  all  within  200m  of  one  another.  The  two 
ADCP  sites  (Al  and  A2)  were  about  29  km  apart,  and  the  distance  between  station  A2 
and  its  nearest  CPIES  location  (Cl)  is  about  12  km. 


Data 

Although  ADCP  recoveries  at  Al  and  A2  were  successful,  recorded  data  could  not  be 
extracted  at  Al  during  Legs  I  and  II  due  to  the  failure  of  memory  cards  inside  the  ADCP. 
The  higher  capacity  memory  cards  were  separately  purchased  (different  from  those 
provided  by  the  manufacturer).  Recording  tests  prior  to  the  deployments  were  successful, 
and  it  is  not  clear  what  caused  this  failure. 
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Figures  20  and  21  show  raw  records  of  temperature,  pressure,  heading,  pitch  and  roll  at 
moorings  A1  and  A2.  Despite  the  proximity  of  A2  positions  for  all  Legs  (within  200  m  of 
one  another),  the  mean  recorded  pressure  difference  between  Legs  I,  HI  and  Leg  II  is 
about  45  dbar.  Moreover,  a  tidal  signal  superimposed  on  a  linear  trend  was  obvious 
during  Leg  II,  while  tidal  signals  were  indistinct  during  Legs  I  and  III.  The  tidal  signal 
was  also  seen  in  the  pressure  record  at  A1  for  Leg  III.  The  pressure  records  were  thought 
to  be  erroneous  at  A2  during  Legs  I  and  III,  presumably  due  to  the  pressure  sensor  being 
deployed  about  90  m  beyond  its  limit.  The  ADCP  (Serial  No.  2937)  used  at  A2  during 
Legs  I  and  III  has  a  200  m-rated  housing.  This  was  upgraded  to  a  500  m-rated  housing  for 
the  deployment  at  a  depth  greater  than  200  m.  The  original  pressure  sensor,  however,  was 
rated  for  200  m,  and  was  not  upgraded.  Apart  from  these  problems,  pressure, 
temperature,  and  current  data  records  are  complete.  The  pitch  and  roll  of  ADCP 
transducers  changed  little  at  both  locations  during  the  mooring  periods.  The  heading 
remained  almost  unchanged  for  the  TRBM-ADCP  package  at  Al.  On  the  other  hand, 
heading  variation  was  large  at  mooring  A2,  especially  during  Leg  III.  During  Leg  II,  the 
heading  changed  abruptly  from  about  250°  to  90°  in  December  2003  and  then  gradually 
increased  to  about  160°.  (Due  to  the  limited  vertical  range  of  WHS300,  current  data 
return  was  reduced  below  50%  at  depths  shallower  than  the  30th  bin  (30m)  at  Al  and  the 
14th  bin  (171m)  at  A2  (Tables  9  and  10). 

Figures  22  and  23  show  the  raw  current  data  from  Al  and  A2,  respectively.  The  first  bin 
depth  at  each  mooring,  determined  by  the  pressure  records,  is  10.1  m  and  6.0  m  above  the 
ADCP  mooring  depth  for  8  m  and  4  m  bin  lengths,  respectively.  The  first  bin  depth  at 
each  mooring  and  on  each  leg  is  listed  in  Table  7. 


ADCP  housing 
'WHS  6  buoy) 


Plastic  buoy 


:oustic  release 


Jhain  weight 


Anchor 


Pop-up  buoy 
Ballast  (main  buoy) 

Cement 


ADCP 


Figure  19.  Two  types  of  ADCP  moorings:  (a)  TRBM-ADCP  package  atAl,  (b)  1-type  subsurface  mooring 

line  at  A2. 
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Figure  20.  Records  of  temperature,  pressure,  heading  (FT),  pitch  (P),  and  roll  (R)  at  mooring  Al.  Tick 

marks  indicate  the  beginning  of  the  month. 


30 


Press,  at  A2  (dbai) 


Figure  21.  Records  of  temperature,  pressure,  heading  (H),  pitch  (P),  and  roll  (R)  at  mooring  A2.  Tick 

marks  indicate  the  beginning  of  the  month. 
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Oct  Nov.  Dec.  Jan.  Feb.  Mar.  Apr.  May  Jun.  Jul.  Aug.  Sep.  Oct  Nov. 
2003  2004 


Figure  22.  Raw  ADCP  currents  observed  at  mooring  Al.  Vertical  axis  is  towards  True  North.  Tick  marks 

indicate  the  beginning  of  the  month. 
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Figure  23.  Raw  ADCP  currents  observed  at  mooring  A2.  Vertical  axis  is  towards  True  North.  Tick  marks 

indicate  the  beginning  of  the  month. 
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Table  9.  Mooring  Al.  Basic  statistics  for  cross-shelf  (Ur,  toward  128°)  and  along-shelf  (Vr,  toward  38°) 
components  of  despiked,  unfiltered  currents,  and  data  return  in  %>.  Max.  and  Min.  denote  recorded 
maximum  and  minimum  values  of  current  components.  The  direction  of  vector  currents  is  measured 
clockwise  from  True  North. 


Bin 

Depth 

(m) 

V  elocity 

Mean 

(cm/s) 

STD 

(cm/s) 

Max. 

(cm/s) 

Min 

(cm/s) 

V  ector  Mean 

Data 

Returned  (%) 

Direction  (°) 

Speed  (cm/s) 

29 

34 

Ur 

-as 

25.7 

1202 

-1203 

305 

3Q6 

508 

Vr 

306 

23.4 

1303 

-80.2 

28 

38 

Ur 

-10 

25.4 

1105 

-1208 

301 

3Q2 

714 

Vr 

302 

22,6 

127.9 

-71,6 

27 

42 

Ur 

-12 

25.1 

1041 

-1206 

358 

208 

826 

Vr 

208 

21,7 

1244 

-65,0 

26 

46 

Ur 

-14 

25,0 

1000 

-1242 

303 

201 

917 

Vr 

200 

21.0 

1209 

-50.4 

26 

60 

Ur 

-15 

24,9 

704 

-1005 

300 

286 

959 

Vr 

206 

20.5 

1104 

-52.4 

24 

64 

Ur 

-16 

24.6 

709 

-96.6 

347 

27.9 

909 

Vr 

27.9 

20.0 

1007 

-44.5 

23 

58 

Ur 

-15 

24.2 

705 

-85.7 

348 

27.2 

907 

Vr 

27.1 

19.6 

900 

-42.1 

22 

62 

Ur 

-15 

24.1 

725 

-84,6 

348 

204 

1000 

Vr 

203 

19,5 

90S 

-41,1 

21 

66 

Ur 

-16 

24,0 

741 

-84,2 

340 

206 

1000 

Vr 

205 

19.3 

940 

-38.0 

20 

70 

Ur 

-16 

23.8 

8Q5 

-84.3 

341 

245 

1000 

Vr 

244 

19.4 

905 

-36.9 

19 

74 

Ur 

-18 

23.6 

815 

-83.2 

336 

232 

1000 

Vr 

202 

19.5 

1005 

-56.5 

18 

78 

Ur 

-20 

23.5 

707 

-82.5 

329 

218 

1000 

Vr 

217 

19.5 

1032 

-63.4 

17 

82 

Ur 

-21 

23,5 

828 

-81,3 

321 

201 

1000 

Vr 

200 

19,4 

1032 

-65,6 

16 

86 

Ur 

-21 

23.3 

800 

-80.6 

316 

106 

1000 

Vr 

104 

19,2 

1003 

-71,0 

16 

90 

Ur 

-20 

22.9 

927 

-79.6 

313 

17.2 

939 

Vr 

171 

18.6 

1031 

-75.8 

14 

94 

Ur 

-19 

22.5 

941 

-76.0 

312 

159 

908 

Vr 

108 

17.8 

800 

-65.0 

13 

98 

Ur 

-18 

22.3 

904 

-80.6 

309 

145 

909 

Vr 

144 

17,2 

823 

-71.1 

12 

102 

Ur 

-16 

22,1 

948 

-81,8 

310 

132 

909 

Vr 

101 

16,7 

815 

-69,6 

11 

106 

Ur 

-14 

21,8 

939 

-80,4 

312 

120 

909 

Vr 

119 

16.3 

8Q4 

-69.5 

10 

110 

Ur 

-13 

21.6 

800 

-76.5 

313 

110 

939 

Vr 

109 

15.9 

8Q5 

-65.9 

9 

114 

Ur 

-12 

21.2 

813 

-76.2 

314 

102 

1000 

Vr 

101 

15.5 

773 

-65.2 

8 

118 

Ur 

-11 

20,9 

736 

-82,0 

312 

04 

1000 

Vr 

04 

15.1 

701 

-62,4 

7 

122 

Ur 

-11 

20,7 

737 

-81,8 

306 

06 

1000 

Vr 

05 

14.6 

720 

-60.8 

6 

126 

Ur 

-13 

20,7 

723 

-82,4 

207 

78 

1000 

Vr 

77 

14.2 

625 

-56.6 

6 

130 

Ur 

-15 

20.7 

749 

-78.7 

206 

71 

1000 

Vr 

09 

13.8 

405 

-50.7 

4 

134 

Ur 

-19 

20.8 

745 

-77.1 

108 

02 

1000 

Vr 

09 

13.5 

400 

-57.0 

3 

138 

Ur 

-24 

20,7 

701 

-74,7 

105 

52 

1000 

Vr 

46 

13,2 

448 

-58,2 

2 

142 

Ur 

-27 

20,1 

731 

-70,3 

-19 

41 

1000 

Vr 

02 

12,6 

4Q6 

-51,0 

1 
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Ur 

-27 

18.5 

610 

-66.0 

-225 

31 

1000 

Vr 

15 

11.0 

375 

-36.1 

34 


Table  10.  Mooring  A2.  Basic  statistics  for  cross-shelf  (Ur,  toward  128  °)  and  along-shelf  (Vr,  toward  38°) 
components  of  despiked,  unfdtered  currents,  and  data  return  in  %.  Max.  and  Min.  denote  recorded 
maximum  and  minimum  values  of  current  components.  The  direction  of  vector  currents  is  measured 
clockwise  from  True  North. 


Bin 

Depth 

(m) 

Velocity 

Mean 

(cm/s) 

STD 

(cm/s) 

Max. 

(cm/s) 

Min. 

(cm/s) 

Vector  Mean 

Data 

Returned  (%) 

Direction 

(°) 

Speed 

(cm/s) 

13 

179 

Ur 

0.5 

15.8 

64.4 

-60.9 

38.7 

46.2 

61 .1 

Vr 

46.2 

25.9 

158.8 

-50.2 

12 

187 

Ur 

-0.3 

15.7 

62.8 

-67.9 

37.6 

44.6 

81 .4 

Vr 

44.6 

26.4 

154.8 

-67.4 

11 

195 

Ur 

-1.0 

15.6 

60.0 

-59.4 

36.7 

42.0 

94.7 

Vr 

42.0 

26.9 

149.9 

-52.2 

10 

203 

Ur 

-1.4 

15.6 

59.3 

-58.4 

35.9 

38.0 

98.9 

Vr 

38.0 

27.8 

144.8 

-52.2 

9 

211 

Ur 

-1.5 

15.6 

69.2 

-55.6 

35.4 

32.6 

99.7 

Vr 

32.5 

28.5 

143.9 

-51.1 

8 

219 

Ur 

-1.5 

15.4 

50.5 

-60.2 

34.7 

26.3 

99.7 

Vr 

26.3 

28.8 

137.6 

-49.9 

7 

227 

Ur 

-1.3 

15.2 

51.1 

-59.3 

34.2 

20.1 

99.7 

Vr 

20.0 

28.3 

131.1 

-47.4 

6 

235 

Ur 

-0.8 

15.1 

54.9 

-52.2 

34.8 

14.4 

99.7 

Vr 

14.4 

27.1 

128.2 

-53.2 

5 

243 

Ur 

-0.2 

15.0 

57.1 

-56.7 

36.5 

9.5 

99.7 

Vr 

9.5 

25.3 

122.6 

-54.6 

4 

251 

Ur 

0.3 

15.1 

63.6 

-60.2 

41.0 

5.5 

99.7 

Vr 

5.5 

23.2 

112.0 

-56.7 

3 

259 

Ur 

0.7 

15.2 

73.0 

-51.5 

54.5 

2.5 

99.7 

Vr 

2.4 

20.9 

98.4 

-56.2 

2 

267 

Ur 

1.2 

15.3 

77.7 

-48.4 

108.4 

1  .3 

99.7 

Vr 

0.4 

18.4 

85.9 

-51.2 

1 

275 

Ur 

1.6 

15.1 

68.6 

-50.1 

137.7 

1  .7 

99.7 

Vr 

-0.3 

15.3 

69.2 

-44.2 
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Data  Processing 


Pressure  and  temperature 

The  accuracy  and  precision  of  the  RDI  temperature  sensor  are  ±0.4°C  and  the  resolution 

is  0.0 1°C.  The  accuracy  and  precision  of  the  RDI  pressure  sensor  are  0.25%  of  the  rated 
sensor  range  and  the  resolution  is  0.01%  of  the  range  (200  m  for  the  sensors  on  ADCP 
Serial  No.  2937  and  2279,  and  500  m  for  sensor  on  ADCP  Serial  No.  3569).  No 
information  is  available  on  the  stability  of  the  sensors. 

Tidal  signals  are  obvious  in  pressure  records  from  both  sites  (except  Legs  I  and  III  at  A2), 
and  there  exists  a  pressure  drift  during  Leg  II  at  A2  while  the  drift  is  indistinct  at  A1 
(Figures  20  and  21).  The  drift  at  A2,  determined  by  fitting  a  linear  curve,  was  0.0028 
dbar/day.  The  record-length  mean  pressures  at  A1  and  A2  are  152.4  dbar  and  286.3  dbar, 
respectively.  Each  ADCP  record  was  used  to  determine  the  depth  of  the  pressure 
transducer  by  converting  the  pressure  record  to  depth  using  a  pre-set  salinity  value  of 
34.0  and  the  following  equation  provided  by  the  manufacturer: 

depth  (in  dm)  =  recorded  pressure  (in  kPa)  x  (1.02  -  [0.00069  x  S]), 

where  S  is  salinity.  The  resulting  mean  depth  of  each  mooring  is  151.7  m  and  285.0  m  at 
A1  and  A2,  respectively.  Since  the  three  deployment  locations  at  A2  are  all  within  200  m 
of  one  another,  representative  depths  of  ADCP  moorings  for  the  entire  observation 
periods  are  taken  to  be  152  m  at  A 1  and  285  m  at  A2,  based  on  the  pressure  records 
during  Leg  III  at  A1  and  Leg  II  at  A2. 

The  procedures  for  despiking  and  low-pass  filtering  of  temperature  records  are  the  same 
as  those  for  the  current  data. 


Current  data 

Data  from  A2  are  combined  for  each  leg  and  for  each  8  m  bin  to  make  a  13-month 
velocity  time  series  at  each  bin  depth.  Gaps  in  the  time  series  at  mooring  A2  occur  for 
two  reasons.  First,  shorter  gaps  exist  because  of  the  time  between  mooring  recovery  and 
re-deployment,  25.5  hours  between  Legs  I  and  II,  and  3.5  hours  between  Legs  II  and  III. 
Second,  longer  gaps  in  the  upper  bin  depths  are  due  to  the  limited  vertical  range  of  the 
WHS300  ADCP.  In  processing  the  current  data,  the  depth  bins  where  the  data  return  is 
less  than  50%  of  the  complete  record  are  excluded.  The  resulting  current  records  are  from 
146  m  to  34  m  depth  at  4  m  intervals  for  Al,  and  from  275  m  to  179  m  depth  at  8  m 
intervals  for  A2.  The  longest  gap  is  77  hours  at  Al,  and  38  hours  at  A2. 

After  the  correction  for  local  magnetic  variation  of  5°W  (http://www.ngdc.noaa.gov/ 
seg/geomag/jsp/IGRFGrid.jsp),  data  spikes  in  U  (eastward)  and  V  (northward) 
components  are  removed  first  by  calculating  the  mean  and  standard  deviation  for  the 
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entire  record  and  eliminating  values  departing  from  the  mean  by  more  than  5  standard 
deviations,  and  then  by  calculating  the  mean  for  every  50  ensembles  (each  ensemble  is  30 
minutes  long,  see  Table  7)  and  eliminating  values  departing  from  this  mean  by  more  than 
4  standard  deviations.  De-spiked  gappy  data  are  then  decomposed  into  cross-shelf  {Ur, 
towards  128°)  and  along-shelf  {Vr,  towards  38°)  components  and  then  de-spiked  again. 
Tables  9  and  10  give  the  basic  statistics  for  Ur  and  Vr  at  each  bin  depth.  In  general,  the 
velocity  decomposition  could  be  based  on  either  the  local  topography  or  the  mean  current 
direction  or  the  direction  of  principal  axis  of  current  fluctuations.  Here,  the  along-shelf 
direction  is  chosen  to  be  perpendicular  to  the  line  of  the  CPIES  array  keeping  in  mind  the 
eventual  calculation  of  volume  transport  with  data  from  both  the  ADCP  and  CPIES 
moorings.  Mean  current  directions  at  the  selected  top  bins  for  A1  and  A2  are  36.5°  and 
38.7°,  respectively,  almost  the  same  as  the  along-shelf  direction  (38°). 

Linear  interpolations  are  applied  individually  to  the  de-spiked  gappy  Ur  and  Vr  time 
series  to  fill  the  gaps.  The  velocity  and  temperature  time  series  are  then  low-pass  filtered 
with  a  2nd  order  Butterworth  filter,  run  forward  and  backward,  with  cutoff  period  of  48 
hours,  and  subsampled  at  12-hour  intervals. 

Figures  24  and  25  show  the  time  series  of  low-pass  filtered  current  vectors  sub-sampled 
every  12  hours  at  moorings  A2  and  Al,  and  Figures  26  and  27  show  low-pass  filtered  Ur, 
Vr,  and  temperature  at  those  two  moorings. 
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Depth  (m) 
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Figure  24.  Vector  time  series  of  low-pass  filtered  12-hourly  currents  at  mooring  Al.  Vertical  axis  is 
towards  True  North.  Tick  marks  indicate  the  beginning  of  the  month. 
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Nov. 


Figure  25.  Vector  time  series  of  low-pass  filtered  12-hourly  currents  at  mooring  A2.  Vertical  axis  is 
towards  True  North.  Tick  marks  indicate  the  beginning  of  the  month. 
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Figure  26.  Mooring  Al.  Time  series  of  along-shelf  (Vr)  (top  panel)  and  cross-shelf  (Ur)  ( bottom  panel) 
components  of  low-pass  filtered  12-hourly  currents,  and  low-pass  filtered  12-hourly  temperature  (middle 
panel).  Tick  marks  indicate  the  beginning  of  the  month. 
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Figure  27.  Mooring  A2.  Time  series  of  along-shelf  (Vr)  (top  panel)  and  cross-shelf  (Ur)  (bottom  panel) 
components  of  low-pass  filtered  12-hourly  currents,  and  low-pass  filtered  12-  hourly  temperature  (middle 
panel).  Tick  marks  indicate  the  beginning  of  the  month. 
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CTD  data 


CTD  data  were  collected  along  the  section  shown  in  Figure  18  during  each  of  four  cruises 
using  a  SeaBird  SBE  911  Plus  CTD  system.  Table  1 1  lists  the  locations  and  depths  of  the 
CTD  casts.  The  CTD  casts  were  taken  at  4  stations  in  October  2003  (during  the  initial 
deployment  cruise)  and  at  1 1  stations  in  May  and  November  2004.  The  locations  of  CTD 
stations  A1  and  A2  and  Cl  through  C6  shown  in  Figure  18  almost  coincide  with  the 
locations  of  the  ADCP  and  CPIES  moorings.  CTD  data  were  taken  at  9  stations  along  the 
PN-line  in  November  2003  on  board  T/V Kagoshima-maru  (Figure  28),  and  the  stations 
are  listed  in  Table  11.  Downcast  CTD  data  were  processed,  binned  into  lm  depth  bins, 
with  parameters  derived  using  the  Seabird  SEASOFT  software  package.  The  CTD  data 
taken  in  November  2003  were  processed  at  Kagoshima  University.  Figure  29  shows  the 
resulting  vertical  sections  of  temperature  and  salinity  with  Figure  29b  being  taken  from  a 
cruise  report  prepared  by  Dr.  Hiroshi  Ichikawa  during  the  November  2003  cruise. 


41 


Table  11.  Locations  of  CTD  stations,  water  depths,  and  CTD  lowering  depths.  The  listed  time,  position  and 
water  depths  are  the  averages  between  the  cast  start  and  end  values. 


Station 

(CTD) 

Date 

Time 

(GMT) 

Position 

Water 

depth  (m) 

CTD 

maximum 

depth  (m) 

Long.  (°E) 

Lat.  (°N) 

A1 

2003/1 0/D  1 

17:31 

126.6856 

28.5581 

154 

153 

58.3 

A2 

M 

15:23 

126.9083 

28.4080 

280 

278 

30.9 

Cl 

n 

14:18 

126.9926 

28.3619 

437 

431 

21.3 

C2 

" 

12:16 

127.1614 

28.2416 

973 

971 

0.0 

FN9 

2003/11/08 

07:00 

126.4617 

28.7635 

124 

145.9 

FN8 

II 

08:22 

126.6130 

28.6644 

139 

127.5 

PN7 

II 

09:36 

126.7620 

28.5701 

190 

109.5 

PN6 

II 

11:04 

126.9185 

28.4772 

273 

91.1 

PN5 

" 

12:43 

127.0720 

28.3818 

842 

72.8 

PN4 

n 

15:08 

127.2204 

28.2901 

981 

55.0 

FN3 

n 

17:23 

127.3692 

28.1890 

1073 

36.6 

PN2 

II 

19:46 

127.5180 

28.0882 

1240 

18.3 

PN1 

" 

21:49 

127.6644 

27.9868 

1303 

0.0 

HI 

2004/05/14 

07:56 

125.9493 

29.0108 

115 

110 

274.7 

H2 

II 

09:50 

126.1884 

28.8501 

104 

96 

245.4 

H3 

2004/05/15 

11:48 

126.4377 

28.7006 

125 

122 

216.0 

A1 

II 

13:40 

126.6708 

28.5563 

147 

144 

188.2 

A2 

" 

17:11 

126.9205 

28.4134 

270 

266 

159.2 

Cl 

II 

18:53 

126.9866 

28.3544 

395 

389 

150.0 

C2 

n 

20:44 

127.1531 

mwcinwi 

00 

’rT 

O, 

951 

128.7 

C3 

II 

13:01 

127.3267 

28.1250 

943 

912 

108.1 

C4 

" 

10:16 

127.5706 

27.9764 

1274 

1264 

79.1 

C5 

" 

07:04 

127.8550 

27.7601 

1073 

1044 

42.3 

C  6 

2004/05/13 

04:14 

128.1989 

27.5321 

833 

821 

0.0 

HI 

2004/11/02 

07:16 

125.9511 

29.0103 

110 

105 

275.1 

H2 

II 

05:42 

126.1882 

28.8485 

103 

99 

245.9 

H3 

2004/11/03 

04:09 

126.4390 

28.6993 

123 

121 

216.4 

A1 

n 

01:55 

126.6713 

28.5556 

144 

142 

188.4 

A2 

n 

14:13 

mw.wKm 

28.4082 

236 

252 

159.9 

Cl 

II 

10:26 

126.9835 

28.3517 

374 

353 

150.7 

C2 

" 

12:02 

127.1532 

28.2338 

929 

896 

129.5 

C3 

2004/11/04 

05:24 

127.3217 

28.1225 

958 

895 

108.9 

C4 

II 

02:57 

127.5725 

27.9699 

1238 

1194 

79.1 

C5 

n 

23:58 

127.8609 

27.7567 

1042 

994 

42.2 

C  6 

" 

21:07 

128.2026 

27.5280 

829 

795 

0.0 
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Figure  28.  Positions  of  CTD  stations  (solid  circles)  occupied  in  November,  2003.  ADCP  mooring  positions 
(Al,  A2)  are  also  shown  (solid  squares).  Bathymetry  is  based  on  ET0P02  in  meter  units. 
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Figure  29.  Temperature  (left  panels)  and  salinity  (right  panels)  along  CTD  sections  measured  four  times  in 
October,  November  2003,  and  May,  November  2004  (from  top  to  bottom). 
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